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ABSTRACT 

We study the time dependent spectra produced via the bulk Compton process by a 
cold, relativistic shell of plasma moving (and accelerating) along the jet of a blazar, 
scattering on external photons emitted by the accretion disc and reprocessed in the 
broad line region. Bulk Comptonization of disc photons is shown to yield a spectral 
component contributing in the far UV band, and would then be currently unobserv- 
ablc. On the contrary, the bulk Comptonization of broad line photons may yield a 
significant feature in the soft X-ray band. Such a feature is time-dependent and tran- 
sient, and dominates over the non thermal continuum only when: a) the dissipation 
occurs close to, but within, the broad line region; b) other competing processes, like 
the synchrotron self-Compton emission, yield a negligible flux in the X-ray band. 
The presence of a bulk Compton component may account for the X-ray properties of 
high redshift blazars that show a flattening (and possibly a hump) in the soft X-rays, 
previously interpreted as due to intrinsic absorption. We discuss why the conditions 
leading to a detectable bulk Compton feature might be met only occasionally in high 
redshift blazars, concluding that the absence of such a feature in the spectra of most 
blazars should not be taken as evidence against matter-dominated relativistic jets. 
The detection of such a component carries key information on the bulk Lorentz factor 
and kinetic energy associated to (cold) leptons. 

Key words: radiation mechanisms: non thermal — scattering — quasars: general — 
quasars: individual: GB B1428-I-4217 — X-rays: general 



1 INTRODUCTION 

High redshift blazars are the most powerful observed radio- 
loud AGN, and among the best candidates to be detected by 
the forthcoming GLAST 7-ray satellite. Recently, one previ- 
ously unidentified EGRET source has been associated with 
the blazar Q0906-693 at a redshift 5.47 (Romani et al. 2004). 
Its spectral energy distribution (SED) resembles the predic- 
tions of the proposed 'blazar sequence' (Fossati et al. 1998; 
Ghisellini et al. 1998), namely the presence of two broad 
peaks in the mm-far IR and the MeV-GeV bands, with 
the high energy component strongly dominating the energy 
output. Furthermore, there are now half a dozen blazars at 
z > i, yet undetected in the MeV-GeV band, but observed 
in the X-ray band with large area instruments (i.e. XMM- 
Newton, see Yuan et al. 2006). All of them exhibit a flat 
(i.e. rising in vF^) spectrum, indicating the energetic domi- 
nance of a spectral component peaking at high frequencies. 
The inferred apparent luminosities reach > 10''^ erg s~^, yet 
the intrinsic ones (i.e. corrected for the effect of relativistic 
beaming) are likely to be only a minor fraction of the jet 
power which energizes the radio lobes. 



Despite the theoretical advances following the discovery 
of the radiatively dominant 7-ray emission in blazars, and 
more than thirty years after the discovery of superluminal 
motion in 3C 273 (Whitney et al. 1971; Cohen et al. 1971) we 
are still unable to pin down the acceleration and coUimation 
mechanism(s) and establish whether the energy carried in 
jets is mainly in the form of Poynting flux or kinetic energy of 
matter and how this might depend on the distance from the 
powering source (see the arguments presented by Blandford 
2002 and Lyutikov 2003 in favour of electromagnetically- 
dominated jets, and Sikora et al. 2005 for a discussion of 
matter vs magnetic jets). 

One of the proposed diagnostics on the issue of jet com- 
position (e.g. Sikora & Madejski 2000) is based on the spec- 
tral signature expected from a matter-dominated jet, in the 
form of an excess of emission at around 1 keV, due to the 
bulk Comptonization by cold leptons of the radiation fields 
produced by the accretion disc and/or reprocessed in the 
broad line region (BLR). This process, first suggested by 
Begelman & Sikora (1987), was then considered by Sikora, 
Begelman & Rees (1994) and Sikora et al. (1997). Moderski 
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et al. (2004) proposed that such a process can generate a 
soft X-ray precursor in the hght curve of blazars. 

No evidence for bulk Comptonization features have 
been detected so far, leading to upper limits on the jet mat- 
ter content (sec above references) . However, there is increas- 
ing recent evidence that the X-ray spectrum of some high 
redshift blazars show a flattening towards low X-ray ener- 
gies and possibly a hump. The favoured interpretation so fax 
has been the presence of intrinsic absorption by warm ma- 
terial, but alternative explanations are viable (sec Fabian ct 
al. 2001b, Worsley et al. 2004). In particular it is possible 
that for the first time we are observing the signature of the 
bulk Comptonization process. 

This motivates the present work which concentrates on 
the (time dependent) detailed shape of the bulk Compton 
spectrum to assess this hypothesis against available and fu- 
ture data, under reasonable assumptions about the pres- 
ence of ambient seed photons and the bulk acceleration of 
the plasma in the jet. Future X-ray missions could pro- 
vide variability as well as polarimetric information in the 
X-ray band, the latter being a further valuable tool to sin- 
gle out emission via bulk Comptonization among competing 
radiation processes (Begelman & Sikora 1987), such as syn- 
chrotron self-Compton (SSC) (Celotti & Matt 1994), and 
external Compton (Poutanen 1994). 



2 MODEL ASSUMPTIONS 

We consider a shell comprising N cold (i.e. non relativistic) 
leptons moving along a jet, accelerating up to a saturation 
speed. During propagation, the leptons scatter ambient pho- 
tons, originating both in a (standard) radiatively efficient 
accretion disc and reprocessed in the gas clouds forming the 
BLR, as expected in powerful quasars. In the following we 
detail the assumptions made in this scenario. 



2.1 Plasma dynamics 

The plasma dynamics, and in particular the dependence 

of the jet velocity on distance z from the power source, 
is not known. Modeling of magnetically accelerated flows 
leads to relatively slow acceleration (e.g. Begelman & Li 
1994), before the flows reaches an asymptotic speed. Here, 
we parametrize the acceleration in terms of the bulk Lorentz 
factor r increasing with distance (i.e. the height above the 
disc) as: 



up to a maximum value Fmax, beyond which F = Fmax. The 
initial bulk Lorentz factor, formally a free parameter, will 
be considered to be 1, unless otherwise specified. The initial 
height Z() (corresponding to the typical distance where accel- 
eration sets in) will be taken of the same order as the inner 
accretion radius ro- To simplify the notation, we hereafter 
set F = F(2) and (3 = f5{z) = (F^ - ly/^/T. 



2.2 Scattering of radiation from the accretion disc 

To model the radiation produced by the accretion disc, we 
simply consider the surface (multi-colour blackbody) tem- 
perature profile for a Shakura & Sunyaev (1973) disc: 
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where r is the radial disc coordinate, rg = GMbh/c^ the 
gravitational radius, Mbh the central black hole mass and 
M is the mass accretion rate. 

At a given height z, the bolometric radiation energy 
density produced by a 'ring' of the accretion disc, as seen in 
the lab frame, corresponds to 
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B{Td) = (JsbTd/ti" represents the frequency-integrated 
black body intensity, ctsb is the Stefan-Boltzmann constant, 
and cos^ = [1 + ^ being the angle between the 

position z along the jet and the emitting ring of the disc. 
In the comoving (primed) frame of the shell this radiation 
energy density is given by: 

^=F^(l-/Jcos.)^^. (4) 
dr dr 

We assume that the scattering process is isotropic, i.e. we 

neglect the angle dependence of the differential Thomson 
cross section, by setting da/dQ. = (7T/(47r). Thus, for an 
optically thin shell (see Sect. 2.5), each electron scatters a 
fraction (tt of this radiation energy density, re-isotropizing 
it in the comoving frame. For N free and cold electrons the 
scattered radiation (at a given z) corresponds to a bolomet- 
ric luminosity (in the comoving frame) 
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dr dr 
An observer located at an angle 6y with respect to the jet 
axis will receive a power 
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defines the relativistic 



where S = r-'^{l - Pcosey) 
Doppler factor. In terms of the spectrum, the observed ra- 
diation will still have a blackbody frequency distribution, 
corresponding to a transformed temperature: 

1 -/3 cos 6* 



TD,ohs{r,z) = To{r) 



1 — /? cos 6\ 



(7) 



Given the total luminosity (eq. 6) and spectrum (eq. 7) , 
the observed spectrum can be normalized by setting 
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(9) 



z2-Hr2 [F(l -/3cos6')]2' 

The observed power emitted at each z can be simply ob- 
tained by integrating over the disc radii r, between ro and 
an outer disc radius (sec Sec. 2.5). 

In the following we will also consider the time integrated 
spectrum, as measured by the observer over an (integration) 
time tohs, given by: 
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2.3 Scattering of broad line photons 

Besides the photons directly originating in the disc, a signifi- 
cant contribution to the soft photon field can come from disc 
photons reprocessed by the gas permeating the BLR (Sikora, 
Begelman & Rees 1994). We account for this contribution as- 
suming that the photon energy distribution follows a black- 
body spectrum peaking at the frequency of the Lyman-a 
hydrogen line, VLa = 2.47 X 10^ ' Hz. This matches the shape 
of this external radiation component in a restricted energy 
range around the Lyman-a line, as seen in the comoving 
fi^ame. 

In fact for BLR clouds distributed in two semi-spherical 
shells (one for each side of the accretion disc) , any monochro- 
matic line is seen, in the comoving frame, within a narrow 
cone of semi-aperture angle l/F along the jet velocity di- 
rection. Accordingly, such photons are blueshifted by a fac- 
tor ranging from T (photons from the border of the cone) 
to 2T (photons head-on). In this (admittedly narrow) fre- 
quency range a monochromatic line transforms into a spec- 
trum oc Vohs^ ■ Although the Lyman-a line represents the 
most important contribution of the BLR seed photons, the 
entire spectrum produced by the BLR is more complex (see 
Tavecchio et al., in prep.). Here, we stress that the blackbody 
assumption represents well the more complex SED from the 
BLR. 

Thus we set an equivalent temperature of the BLR spec- 
trum: 



Tblr = ^ - 4.23 X 10^ K. 



(12) 



In the shell comoving frame this temperature is seen 
blueshifted by a factor ~ F, while, after scattering, the trans- 
formation in the observer frame introduces an additional 
blueshift, corresponding to a factor 5: therefore TBLR.obs ~ 
F^Tblr. 

The bolometric luminosity resulting from the scattering 
of BLR photons, at a given distance z, is given by 



aTcAf^7BLR^25^ 



(13) 



where ?7blr is the energy density of BLR photons, which is 
considered constant within the typical distance of the BLR, 
defined by a radius -Rblr- For a luminosity I/blr of broad 
line photons, the observed spectrum will have a blackbody 
shape, corresponding to the temperature TBLR.obs: 



2h 1^ 

exp (/li^/fcTBLR,obs) - 1 



(14) 



'-BLR-'BLR^ 



which corresponds to the integrated luminosity given by 
Eq. 13. 



2.4 Disc vs broad line photons 

Since the highest temperature is reached in the innermost 
parts of the accretion disc, the angle between the energetic 
disc photons and the moving shell is typically small. It is 
relatively large only in the vicinity of the disc, but there the 
shell has not yet reached high bulk Lorentz factors. As a 
consequence (see Eq. 7), the observed typical frequency of 
the scattered disc photons will not increase much (for ^ ~ 
and 9^ ^ (3 such photons are even redshifted). 

On the contrary, photons from the BLR are always seen 
head-on, and as such are maximally blueshifted (i.e. by the 
factor V5). For this reason (as shown below) it will be this 
bulk Compton component which can provide a relevant con- 
tribution in the X-ray spectrum of powerful blazars. 



2.5 Model parameters 

In the following we quantitatively specify the choice of model 
parameters adopted in the spectral calculations presented in 

Section 3. 

Disc luminosity — For powerful blazars we assume a typ- 
ical black hole mass Mbh = lO^MgM© and a disc accretion 
rate close to the Eddington one. Therefore 



-/^disc ~ -T/Edd ~ lO'^^Mg erg s \ 



(15) 



and the typical sizes of the inner accretion disc and jet are 
assumed zo = ro = 6rg. The outer accretion disc radius is 
set to lO"* rg. 

Broad line radiation — We model the broad line region 
as a thin spherical shell at distance JJblr- The disc radia- 
tion reprocessed by clouds in the BLR is assumed to have a 

bolometric luminosity Lblr = O.lLdisc- 

It should be noted that the BLR may not be the only 
contributor to seed photons, besides those directly from the 
accretion disc. A molecular torus (emitting infrared pho- 
tons; e.g. Blazejowski et al. 2004), and intra-cloud plasma 
(scattering and isotropizing both disc and also beamed syn- 
chrotron jet radiation; e.g. Ghisellini & Madau, 1996) may 
also provide soft photons. However we neglect such compo- 
nents here, as their contribution is very uncertain. 

Number of electrons — We consider a shell of plasma, 
of some initial size equal to the inner size of the accretion 
disc, ro. As first order estimate on how many particles can be 
present in such a shell, we consider the number of relativistic 
electrons needed to originate the radiation dominating the 
SED (i.e. produced in the jot regions whore most of the 
dissipation occurs) and, assuming that the total number of 
particles is conserved, we treat it as a lower limit to the total 
number of electrons present in the shell. 

This limit is clearly model dependent. We will adopt 
the most common interpretation that the high energy radi- 
ation in the spectrum of powerful blazars is produced via 
inverse Compton emission of a relativistic distribution of 
leptons on the external photon fields discussed in Sec. 2.2 
and 2.3. In this scenario the inferred number density of rel- 
ativistic particles relies on estimates of the magnetic field 
and seed photons energy density required to give raise to 
the observed luminosity. On the other hand, these quanti- 
ties are quite tightly constrained by fitting the whole SED 
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- both the synchrotron and inverse Compton components. 
Furthermore for powerful blazars, the electron distribution 
is basically in the complete cooling regime (namely, elec- 
trons of all energies cool within a light crossing time) , which 
constrains the low energy particle end. 

For a reference value of A'^ = 10^'*Ai'54 particles (see Sect. 
3), which initially occupy a volume of typical size ro — /2rg 
(2rg being the Schwarzschild radius) , the corresponding ini- 
tial optical depth is ro = aTro{N/ro) ~ 7 AN54 / if Mgf . At 
the base of the jet, the shell can be optically thick for scatter- 
ing (depending on /): in the calculations we rescale the bulk 
Compton emission by the factor max(l,T), which accounts 
for the reduction in the scattered flux when r > 1 (self- 
shadowing). However, for a conical jet the region rapidly 
becomes transparent as r oc independently of the shell 
width. After transparency is reached all electrons contribute 
to the emission and the process becomes largely independent 
of the specific geometry (shell or spherical) , and the expres- 
sions given in Sect. 2.4 are appropriate. 

Plasma dynamics — As mentioned above we parametrize 
the uncertain dynamics as r(2:) = (z/zo)". In the following 
examples we adopt a = 1/2, which, for Fmax ~ 10, corre- 
sponds to reach the asymptotic bulk Lorentz factor at dis- 
tances z = zoF^ax ~ Bi few hundreds r-g. In the context of 
the internal shock scenario (Rees 1978; Sikora, Begelman & 
Rees 1994; Ghisellini 1999; Spada et al. 2001) this is approx- 
imately the distance where shell-shell collisions take place, 
i.e. where most of the dissipation occurs. For a = 1 (i.e. more 
rapid acceleration) the results are qualitatively the same, 
but the observed time evolution is quicker. 



3 RESULTS 

Fig. 1 shows a sequence of bulk Compton spectra resulting 
from scattering disc and line photons. The sequence refers 
to spectra at different times (in the observer frame) from 
the beginning of the plasma acceleration at zo = 6rg. Each 
spectrum corresponds to the average luminosity received in 
3 hr of integration time, toba, again as measured in the rest 
frame of the source. The acceleration (F oc z^^^) terminates 
when Fmax = 15 is reached. For comparison the BLR size is 
-Rblr = 10^* cm ~ lO^'rg. 

The total number of cold particles is set at ~ 10^^, 
and the observer is located at a viewing angle 9^ = 3° . 

The figure shows separately the contribution from scat- 
tering of disc (continuous lines) and BLR (dashed lines) pho- 
tons. For clarity the top panel shows the evolution of the 
spectrum in the rising phase of the disc-scattered compo- 
nent (roughly until the Lorentz factor reaches Fmax), while 
the bottom panel shows its decline. 

Such evolution corresponds to the following phases: 

• At early times, the bulk Lorentz factor is small, with 
a correspondingly low scattered luminosity, from both disc 
and line photons. The former component is more prominent 
as the radiation energy density of the disc photons is larger 
at small z. The different slope of the disc-scattered radiation 
with respect to the multi-color black body is due to the 
dilution of photons produced at larger disc radii, weakly 
contributing to the total energy density. 
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Figure 1. Temporal evolution of the bulk Compton spectrum. 
Long dashed line: accretion disc spectrum. Solid hnes: bulk Comp- 
ton spectra from scattered disc photons. Dashed Unes: bulk 
Compton spectra from scattered broad line photons. Each spec- 
trum corresponds to the average luminosity, assuming 3 hours of 
integration time, produced by the shell as it moves from zq to the 
BLR, assumed to be at -Rblr = 10^* cm. The top panel shows 
the evolution for the first 21 hr (observer time) during which the 
contribution from bulk Comptonization of disc photons increases; 
the bottom panel corresponds to the time interval between 21 and 
51 hr, during which the same component decreases in flux. Notice 
that the luminosity and spectrum from the bulk Comptonization 
of BLR photons remain constant after ~ 21 hr (bottom panel), 
corresponding to the fact that F has reached its maximum value. 

• As the jet accelerates (before reaching Fmax), the BLR- 
scattered component increases more rapidly than the disc- 
scattered one. This is because the radiation energy density of 
the BLR photons is constant for a stationary observer inside 
it, but increases like F^ in the shell comoving frame. The 
radiation energy density produced by the inner disc, instead, 
decreases for a stationary observer, and even more so in the 
shell frame (since the photon directions form small angles 
with respect to the shell velocity). Only photons produced 
at large disc radii are blueshifted, and make up the dominant 
part of the disc radiation energy density. 

• When the shell reaches its maximum Lorentz factor (at 
tobs ^ 21 hr), it is already relatively far from the disc. In 
the comoving frame a larger and larger portion of the disc 
is seen at small angles (i.e. in the comoving frame photons 
are seen redshifted). 

Correspondingly, the disc-scattered radiation intensity 
decreases (bottom panel of Fig. 1), while the BLR line- 
scattered component remains constant, and dominates. 

• As mentioned in Sec. 2.4, the typical energies of the disc 
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Figure 2. The bolometric bulk Compton luminosity from the 
scattering of disc photons (top panel) and the frequency peak of 
this component (bottom panel) as a function of the height of the 
shell above the disc, z. Broad line photons are neglected here. 
Thick (black) lines show the contribution from soft photons pro- 
duced in the entire disc, while thin lines represent, respectively, 
the contribution from 6-60 rg (blue); 60—600 rg (green) and 600— 
6000 rg (red). Note that the three regions of the disc similarly 
contribute to the total luminosity. Solid and dashed lines corre- 
spond to viewing angles 6^ = 0° and 6^ = 3° , respectively. We 
assumed Ldisc = 10*^ erg s~^, A/bh = W^Mq. 



and BLR line scattered components are different. The disc- 
scattered hump peaks in the far UV, in the spectral region 
currently most difficult to observe. This is because energetic 
disc photons reach the shell at small angles, except at very 
early times when however the Lorentz factor is still low. 
Therefore they are not significantly boosted. Line photons, 
instead, always intercept the shell head on (in the comoving 
frame), and thus are maximally boosted. The peak frequency 
of this component depends only on the value of F that the 
shell has at a given z (and of course on the viewing angle). 

To investigate the last item better, which is crucial in 
assessing the detectability of bulk Compton emission, we 
show in Fig. 2 the bolometric luminosity and the peak fre- 
quency of the disc-scattered spectrum as a function of the 
shell-disc distance, z. Solid and dashed lines correspond to 
dv = 0° and 9v = 3°, respectively. The contribution of seed 
photons from the inner (6-60 rg), central (60-600 rg) and 
outer (600-6000 rg) parts of the disc are shown separately: 
the inner part is relevant only when the shell is close to the 
disc, while the other part mostly contributes when the shell 
moves beyond a few (~ 10) rg. The curves change behavior 
when F reaches Fmax (or z > -Rblr), which in this particular 
example occurs at z ~ 2 x lO'^^ cm. 
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Figure 3. The bulk Lorentz factor F, Doppler factor 5, and ob- 
served (rest frame) time to^s (the latter expressed in hours), as 
a function of the location of the shell above the disc. The input 
parameters are the same used in the previous case, except for 
^disc = 4 X 10""" erg s"! (for a black hole of Men = 4 X IO^Mq); 
^BLR = 1.2 X 10^* cm; a total number of (cold) particles 
Af = 8.6 X 10^4 and a viewing angle 0v = 3°. The vertical dashed 
lines corresponds to the start of our calculations (i.e. 6rg) and 
the assumed location of the BLR. 



3.1 The case of GB B1428+4217 

To illustrate the possible relevance of the bulk Compton pro- 
cess for the interpretation of the X-ray spectra of blazars, 
which motivates this work, we consider the case of GB 
B1428-I-4217, a powerful blazar at redshift z = 4.72. This 
is one of the best studied high redshift (z > 4) blazars, the 
spectrum of which has been observed in X-rays by ROSAT 
(Boiler et al. 2000), BeppoSAX (Fabian et al. 2001a) and 
XMM-Newton (Worsley et al. 2004; Worsley et al. 2006). 
These observations showed that the X-ray spectra flat- 
ten towards low energies with respect to a higher energy 
power law. This behaviour has been interpreted as likely 
due to (warm) absorbing material (with columns exceeding 
lO^'^ cm"'^). Alternative interpretations are however possi- 
ble (Fabian et al. 2001b): a flattening of the low energy part 
of an inverse Compton component could be due to a low 
energy cut-off in the electron population and/or a sharply 
peaked soft (seed) photon distribution. The sharpness of the 
X-ray break (occurring over a range of only a few keV in 
the rest-frame) is however difficult to explain within these 
models. 

Here we consider the possibility that the flattening and 
(possibly) a hump in the soft X-ray spectrum (see Worsley 
et al. 2004, Yuan et al. 2005) reveal the presence of emission 
via bulk Comptonization. For this specific example we con- 
sider the parameters which allow to reproduce the high en- 
ergy spectral component on the basis of the external inverse 
Compton process (e.g. in Ghisellini, Celotti & Costamante 
2002). The SED can be satisfactorily fitted as emission of 
a shell with cross section radius ~ 2 x lO^'^ cm and width 
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Figure 4. Illustration of how the bulk Compton process can ex- 
plain some details of the X-ray spectrum of a powerful blazars, 
GB B1428+4217. The X-ray data are from Beppo-SAX (red), 
XMM- Newton (cyan) and ASCA (blue). The power law is as- 
sumed to start from the top of the two Comptonization humps, 
corresponding to assume that a number of leptons similar to those 
'cold' are accelerated to relativistic energies and form a power law 
distribution extending down to 7min = 1- Top panel: spectrum 
predicted during 3 hours of exposure time starting 100 hours af- 
ter the beginning of the acceleration. Bottom panel: the same, 
but after 110 hours. Due to the hardness of the power law, the 
contribution to the total power law of disc and broad line pho- 
tons is almost equal even after 100 hours, despite the fact that 
the bulk Compton hump from disc photons is less powerful than 
the one from the broad line photons. All frequencies are in the 
rest frame of the source. 

(as measured in the comoving frame) AR' = 1.2 x 10^^ cm, 
moving with a bulk Lorentz factor F — 17. The inferred mag- 
netic field is _B = 11 G and the (comoving) number density 
of relativistic particles ~ 5.9 x 10^ cm~^, corresponding to a 
total number of relativistic particles A'^rci ~ 8.8 x 10^^. The 
observer is at 6^ = 3° . 

To account also for the soft X-ray flattening and hump, 
we assume that again the plasma accelerates from zo = 6rg 
according to F = (z/zq)^^^, up to Fmax = 17. The dy- 
namics of the shell with respect to the observer time corre- 
sponds to the case presented in Fig. 3. The disc luminosity 
is Ldisc = 4 X 10*^ erg s~\ for a Mbh = 4 x 10® Mq. The 
BLR, located at i?BLR ~ 1-2 x 10^* cm, reprocesses 10 per 
cent of the disc radiation. The resulting spectra are shown 
in Fig. 4, where we report different X-ray data sets, and 
optimize the parameters to reproduce the BeppoSAX (top 
panel) and 'KMM-Newton (bottom panel) states. 

The overall SED clearly depends on the relative nor- 



malisation of the relativistic and cold electron distributions. 
The power law components in Fig. 4 (dotted lines) extend 
from the peaks of the two bulk Compton components: this 
corresponds to a similar number of electrons in the two par- 
ticle distributions 0, namely a total number of (cold) par- 
ticles A'^ — 4.0(1.3) X 10'''^ (for the top and bottom panels, 
respectively). 

The normalizations of the relativistic and cold particle 
distributions have been assumed to be equal, and the two 
populations have been considered independently. Clearly, a 
self-consistent treatment should in principle account for the 
acceleration and cooling processes to determine the relative 
normalization (in number and energy) as a function of time. 
The assumption made here corresponds to about half of the 
particle present at the dissipation site to be accelerated. Al- 
ternatively this could correspond to a situation where all 
leptons are cold for half of the exposure time, and relativis- 
tic for the other half. 

The spectra shown in Fig. 4 refer to the average lumi- 
nosity detected in three hr of integration time (rest frame), 
starting at ti — 100 hr (top panel) and t2 =110 hr (bot- 
tom panel) with respect to the beginning of the acceleration 
(z — Zo). These times correspond to locations of the shell 
quite close to Rblh (see Fig. 3). 

After ti, the luminosities [uL^) of the bulk Compton 
spectra of disc and line-scattered photons are similar, while 
at t2 the contribution from the line photons dominates. 
In fact, by t2 the shell has already reached its maximum 
Lorentz factor, quite far from the inner parts of the disc. 

As already shown (see Fig. 1), the bulk Compton spec- 
trum of the disc photons, at all times, peaks in the far UV, 
and is thus undetectable. However, its luminosity is highly 
time-dependent (see the top panel of Fig. 1) , exceeding the 
bulk Compton luminosity off line photons at times t 
(rest frame) for the specific example shown here. 

As a consequence should dissipation - namely acceler- 
ation of electrons to relativistic energies - occur at lower 
height above the disc, we should have a non-thermal spec- 
trum dominated by the scattering of the disc, not line, pho- 
tons. This would imply that the X-ray excess due to bulk 
Compton off the line photons would become undetectable, 
being exceeded by non-thermal radiation. This can be seen 
by comparing the spectra in the top and bottom panels of 
Fig. 4. Early dissipation, however, faces a severe difficulty 
in accounting for the SED, since the source becomes opaque 
to the 7-7 process. The emission of the created pairs 

would contribute at frequencies between the two broad peaks 
of the SED, resulting in an overproduction of X-rays with 
respect to the observed level (Ghisellini & Madau 1996). 

Note that - as shown in Fig. 4 - the power law contin- 
uum from disc-scattered photons dominates over a limited 
frequency range even if the corresponding bolometric en- 
ergy density is smaller than that of line-scattered photons. 
For non thermal scattered radiation with a hard spectrum 
(a < 1) the typical frequencies of the seed photons are im- 
portant: disc radiation should not be discarded as main con- 



^ For the relativistic particle distribution we assume N('^) oc 
^-(2a+i) extending down to 7inin ~ !> where a is the spectral 
index, L{v) Qcy~°' . 



Bulk Comptonization spectra in blazars 7 



tributor to the seed photon field, on the basis of its bolomet- 
ric energy density only (see Dermer & Schlickeiser 1993). 

In the observationally interesting X-ray range (i.e. 0.1- 
10 keV, that translates into 0.5-50 keV rest frame for sources 
at 2 ~ 4), the spectrum can be approximated as the sum of 
a power law and a blackbody. The latter component would 
correspond to a hump in the mid X-ray band, accompanied 
by a flattening towards soft energies. 

The bulk Compton feature is expected to be time- 
dcpcndcnt, and generally transient. The necessary condi- 
tions for it to be detectable are: i) dissipation must occur far 
from the accretion disc, but within and near to the BLR; ii) 
the shell should move close to its maximum speed; iii) other 
competing processes, such as SSC, should not dominate over 
the bulk Compton and external inverse Compton processes 
(see Sect. 4). 



4 DISCUSSION 

We have shown that bulk Comptonization of the radiation 
from an accretion disc around a supermassive black hole, al- 
though important in a bolomctric sense, would be currently 
unobservable, since its spectrum peaks in the far UV. On the 
contrary, the bulk Compton process on photons produced in 
the broad-line region can give raise to an observable feature 
which can be represented by a blackbody. The entire (bulk 
Compton-l-relativistic Compton) spectral shape, can then be 
approximated as blackbody-fpower law. 

The bulk Compton feature is however difficult to detect, 
even in the most favorable case. Firstly, in general it can be 
a transient feature, if the plasma injection or dissipation in 
the jet is not continuos. Such possibilities are in general sup- 
ported by the strong variability across the electromagnetic 
spectrum of blazars. Above we mimicked these possibilities 
assuming that the dissipation region is a shell moving along 
the jet, thus giving raise to a bulk Compton feature when the 
shell is close to the BLR. However, it is well possible that 
a 'cold' component is in fact more continuous. In such a 
case the bulk Compton component would correspond to the 
'integration' of the spectra shown. Secondly, the radiation 
has to be strongly beamed, namely the jet has to be highly 
relativistic and the angle with the line of sight small. A 
third condition for a bulk Compton feature to be detectable 
is that other emission processes do not mask it. The most 
likely component that could hide it is that produced by SSC. 
It should be noted that the estimated spectrum (see Fig. 4) 
does not include emission via the SSC process, as this would 
require to assume - as another free parameter - the relative 
intensity of the external soft photon and magnetic fields. 
The SSC could in principle provide a significant contribu- 
tion at the UV and softest X-ray energies even in powerful 
sources. However, from the observational point of view, the 
SSC component contributes the least to the X-ray spectrum 
of the most powerful sources, since among blazars they show 
the flattest X-ray spectra and the largest "Compton dom- 
inance" , i.e. the ratio of the high vs low energy component 
increases eis the bolometric luminosity does (Fossati et al. 
1998). This can been interpreted (GhiscUini ct al. 1998) as 
due to the increased importance of the external seed radi- 
ation fleld relative to the locally produced synchrotron one 
with increasing source power. 



This contrasts with the finding that the radiation en- 
ergy density of BLR photons decreases or remains constant 
in higher luminosity objects, according to the relations re- 
ported by Kaspi et al. (2000) and Bentz et al. (2006), respec- 
tively. On the other hand, BLR photons might not be the 
only contributors to the seed field: the BLR clouds them- 
selves and/or some intra-cloud scattering material cam en- 
hance and isotropize beamed synchrotron radiation from the 
jet, the external radiation (Ghisellini & Madau 1996) and 
disc radiation (Sikora, Bcgclman & Roes 1994). The presence 
and the optical depth of such scattering material could de- 
pend on the disc luminosity: more massive outflows or winds 
from the disc could be favoured at accretion rates near the 
Eddington one. The emission lines and thermal optical-UV 
component in high rcdshift blazars indeed indicate accretion 
luminosities close to the Eddington ones even for black hole 
masses exceeding IO^Mq. 

If the Lyman-a hydrogen line provides most of the seed 
photon field, the detection of the bulk Compton signature 
would constitute a powerful diagnostic to estimate the prod- 
uct sr. Another independent estimate of T or 6, would then 
allow to estimate both the jet speed and viewing angle. For 
instance, if the source shows supcrluminal motion with an 
apparent speed /3a (and assuming F does not change) we 
have: 



rs 



/3a 



sin 

- PcosB 



V5 



(16) 



(17) 



where v-p is the peak of the bulk Compton spectrum. The 
above two equations can be solved for 9^ and F: 



(18) 



(19) 



[2i/p/i/L. -1-/3I]'/'' 

The other key diagnostic which can be inferred from a 
bulk Compton feature is the kinetic energy associated to the 
cold scattering material, which is independent of the mag- 
netic field intensity and the filling factor of the scattering 
plasma. 

In the case of GB B1428-I-4217, the above interpretation 
of the SED implies Skin = NVrUpC^ ~ 10(3.4)1.9 x 10^^ erg, 
for a proton-electron plasma. If the width of the shell con- 
taining this mass is of the order of the initial jet size (in the 
observer frame, i.e. Aii ~ 6rg ~ 3.6 x 10^^ cm), then the 
corresponding kinetic luminosity is Lkin = -Ekinc/(6rg) ~ 
8.6(2.8) X 10*^ erg s~^, corresponding to about the lumi- 
nosity radiated by the accretion disc. Lkin is clearly less if 
there are electron-positron pairs; it should be noticed that 
in the case of a pair-dominated plasma, the jet could be 
significantly decelerated because of the drag due to the bulk 
Compton, but even more so when dissipation occurs and 
also the 'rocket' effect would come into play. 

Note also that Lkin oc 1/A7?. AT? is here a free, un- 
known parameter and indeed a different assumption on Aii 
has been made by Ghisellini et al. (2003). However, despite 
this uncertainty on Lkin, the effective L^kin is in fact substan- 
tially independent of such choice. Fitting the overall SED 
would also allow us to estimate the magnetic field intensity 
and hence the power in Foynting flux. To this aim we will 
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examine the X-ray data and SED of powerful, high redshift 
blazars with indications of a bulk Compton feature (Fabian 
et al., in prep.). Interestingly, in agreement with the pre- 
diction of this scenario, indications of soft X-ray flattening 
appears to be connected with the source X-ray luminosity 
(Yuan et al. 2006). 

Future large area X-ray instruments might even pro- 
vide time dependent X-ray spectra from which in principle 
should be possible to constrain the jet acceleration. Nev- 
ertheless, different shells might contribute to the observed 
spectrum and thus the deconvolution from X-ray variability 
to plasma acceleration might be much more complex than 
outlined here. 



5 CONCLUSIONS 

The study of the bulk Compton spectrum produced by cold 
plasma accelerating in a blazar jet showed that: 

1) If dissipation (i.e. acceleration/injection of relativistic lep- 
tons) occurs during the (early) shell acceleration phase, the 
radiation from the disc can provide a relevant source of pho- 
tons for the relativistic external Compton process. Broad 
line radiation becomes competitive and dominant for such 
process at a shell-disc distance comparable to the broad 
line region. The dominance between the two components 
depends on the spectral index of the relativistic Compton 
spectrum. 

2) The bulk Comptonization of the radiation field originat- 
ing in the accretion disc produces a component peaking in 
the far UV, thus most difficult to observe, and its contribu- 
tion to the soft X-ray band is negligible. 

3) Instead the bulk Compton scattering of broad line pho- 
tons may give an observable feature in the soft X-ray band 
of the most powerful quasars, if the dissipation takes place 
close to, but within, the broad line region radius when the 
shell has reached high bulk Lorentz factors (and for small 
viewing angles). 

4) The transient nature of the bulk Compton feature and 
the requirement that other radiation processes (most no- 
tably, the synchrotron self-Compton process) do not dom- 
inate over it, imply that the bulk Compton component is 
rarely detectable. Its absence cannot thus be considered as 
evidence that jets are not matter dominated. 

5) On the other hand, the detection of this feature allows two 
key quantities to be estimated: i) the bulk Lorentz factor, 
and ii) the amount of cold leptons in the shell and in turn 
the corresponding kinetic power. 

6) In principle the time dependent behavior of the bulk 
Compton component would trace the jet acceleration profile. 

7) The bulk Compton contribution would result in a flat- 
tening of the X-ray spectrum towards softer energies, which 
could mimic absorption, and possibly a hump. Indications 
of such features have been found in the spectra of some 
powerful high redshift blazars. An re-analysis of their X-ray 
spectra and SED will be presented in a forthcoming paper. 
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